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We perform a manifestly gauge-independent analysis of the vacuum stability in the standard 
model including two-loop matching, three-loop renormalization group evolution, and pure QCD 
corrections through four loops. All these ingredients are exact, except that light-fermion masses 
are neglected. We in turn apply the criterion of nullifying the Higgs self-coupling and its beta 
function in the modified minimal-subtraction scheme and a recently proposed consistent method for 
determining the true minimum of the effective Higgs potential that also avoids gauge dependence. 

Exploiting our knowledge of the Higgs-boson mass, we derive an upper bound on the pole mass of 
the top quark by requiring that the standard model be stable all the way up to the Planck mass 
scale and conservatively estimate the theoretical uncertainty. This bound is compatible with the 
Monte Carlo mass quoted by the Particle Data Group at the 1.3 ct level. 

PACS numbers: 11.10.Gh, 11.10.Hi, 14.65.Ha, 14.80.Bn 


The standard model (SM) of elementary particle 
physics has been enormously consolidated by the discov¬ 
ery ra at the CERN Large Hadron Collider of a new 
weak neutral resonance that, within the present experi¬ 
mental uncertainty, shares the spin (J), parity (P), and 
charge-conjugation (C) quantum numbers J PC = 0 ++ 
and the coupling strengths with the SM Higgs boson H , 
in the absence of convincing signals of new physics be¬ 
yond the SM. Moreover, its mass of (125.7 ± 0.4) GeV 
[2] falls well inside the Mp range predicted within the 
SM through global analyses of electroweak (EW) preci¬ 
sion data [2]. Besides completing the SM particle multi- 
plet and confirming the Higgs mechanism of mass gener¬ 
ation via the spontaneous breaking of the EW symmetry 
proposed by Englert, Higgs (The Nobel Prize in Physics 
2013), and Brout, this groundbreaking discovery also has 
fundamental cosmological consequences by allowing con¬ 
clusions regarding the fate of the Universe via the anal¬ 
ysis of the vacuum stability [3j. In fact, owing to an 
intriguing conspiracy of the SM particle masses, chances 
are that the Higgs potential develops a second minimum, 
as deep as the one corresponding to the vacuum with ex¬ 
pectation value (VEV) v = 2 _1 / 4 Gj)V 2 = 246 GeV in 
which we live, at a field value of the order of the Planck 
mass Mp = 1.22 x 10 19 GeV [Tj [5j. This would im¬ 
ply that the SM be stable all the way up to the energy 
scale where the unification with gravity is expected to 
take place anyways, which would diminish the necessity 
for grand unified theories at lower scales. EW symme¬ 
try breaking might thus be determined by Planck-scale 
physics [5], and the existence of a relationship between 
Mp and SM parameters might signify a reduction of fun¬ 
damental couplings. Of course, experimental facts that 
the SM fails to explain, such as the smallness of the neu¬ 
trino masses, the strong CP problem, the existence of 
dark matter, and the baryon asymmetry in the Universe, 
would then still call for an extension. 

Obviously, the ultimate answer to the existential ques¬ 


tion whether our vacuum is stable or not crucially de¬ 
pends on the quality of the theoretical analysis as for 
both conceptual rigor and high precision, and it is the 
goal of this Letter to significantly push the state of 
the art by optimally exploiting information that has be¬ 
come available just recently. The technical procedure 
is as follows. The set of running coupling constants, 
including the SU(2)/, U(l)y, and SU(3) C gauge cou¬ 
plings g(n), g'(g ), and g s {^), respectively, the Higgs self¬ 
coupling A(/u), and the Yukawa couplings 2 //(/r), of the 
full SM are evolved in the renormalization scale /i from 
/r thr = 0(v) to g CT1 = 0{Mp) using the renormalization 
group (RG) equations. The beta functions appearing 
therein take a simple polynomial form in the modified 
minimal-subtraction (MS) scheme of dimensional regu¬ 
larization. They are fully known through three loops [6] 
in the approximation of neglecting the Yukawa couplings 
of the first- and second-generation fermions, and the ones 
of g s (.71 and y q [8] also at the four-loop order 0{a ^), the 
latter being given by the quark mass anomalous dimen¬ 
sion. The initial conditions at fj = /z thr are evaluated 
from the relevant constants of nature, including Som- 
merfeld’s fine-structure constant axh defined in Thomson 
scattering—or, alternatively, Fermi’s constant Gf —, the 
strong-coupling constant a[ 5 \Mz) at its reference point 
in QCD with nj = 5 active quark flavors, and the phys¬ 
ical particle masses ( i = W, Z, H , /) defined via the 
propagator poles, taking into account threshold correc¬ 
tions [9], which are fully known through two loops [si ao}- 
fl4] and, for g s and y q , also at 0(ai?) [EE lU] and even at 
O(cts) [EL M- Although self-consistency requires that 
n-loop evolution is combined with (n — l)-loop match¬ 
ing, we, nevertheless, include the additional information 
ELI ag in our default predictions. There are two ap¬ 
proaches to the threshold corrections in the literature 
that differ in the definition of the MS VEV In the 
first one [TOl lllj . v(n) is fixed to be the minimum of the 
effective Higgs potential V e s(H) in the Landau gauge and 
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is thus gauge dependent m- A solution to this problem 
has recently been proposed in Ref. [20]. In the second 
approach El nun, the adjustment of the VEV is only 
done for the bare theory, yielding v° = i/—(m^) 2 /A°, 
with being the mass of the complex scalar doublet 
d>, or, equivalently, 



in terms of basic parameters of the broken phase [9]. The 
linear term in the bare Higgs potential is then quenched 
and cannot serve as a tadpole counterterm, so that the 
tadpole contributions, which carry gauge dependence, 
need to be properly included order by order [3]. Upon 
MS renormalization, taking Eq. ([l) with the superscripts 
0 dropped to be exact, v(g) and all the basic parameters, 
including A (g), are manifestly gauge independent to all 
orders. Consequently, the twofold vacuum stability con¬ 
dition 0, 

A( M cri ) = ^(/x cri ) = 0, (2) 


which fixes a second minimum that is degenerate with the 
first one, has gauge-independent solutions for the critical 
ultrahigh scale g crl and one free basic parameter, which 
we take to be M t crl , the upper bound on the top-quark 
pole mass M t , which is much less precisely known than 
M h 0. For comparisons with the literature, we also de¬ 
termine the Mh lower bound M^ 1 sloppily using as in¬ 
put the mass parameter M t MC [2j that is extracted from 
experimental data using Monte Carlo event generators 
merely equipped with leading-order (LO) hard-scattering 
matrix elements. The results for g crl obtained together 
with M t crl and Mg' 1 are denoted as gf 1 and //A, respec¬ 
tively. While the criticality condition in Eq. pi) carries 
a very simple physical meaning and is straightforward to 
solve numerically, it is slightly scheme dependent. To as¬ 
sess this scheme dependence, we compare the results for 
/r™ and with i = t,H with those obtained applying 
the consistent approach of Ref. |20j . in which V e ff(H) is 
reorganized in powers of h, so that its expansion coeffi¬ 
cients are gauge independent at its extrema [2Tj . Specif¬ 
ically, this amounts to solving 


A = 


1 


256?r 2 
+ 2 g' 4 


(g'+g'y l-31n 


g 2 + g ' 2 


l- 3 ln^_ ) -48 y\ 


1-lnh 


(3) 


which follows from dV^f 1 (g CT1 )/dH = 0, for the mini- 
mum H = g CI1 of V^P(H) and requiring that, at next-to- 
leading orderjNLO), V^° =y^ s °(g cli ) + V^ LO {g cli ) > 
0 for M t < M t crl or Mh > M|p, which is conveniently 
achieved in the Landau gauge pup- 
We adopt the input values G F = 1.1663787(6) x 
1(T 5 GeV" 2 , af\M z ) = 0.1185(6), M w = 

80.385(15) GeV, M z = 91.1876(21) GeV, M H = 


125.7(4) GeV, Af t MC = 173.21(87) GeV, and M b = 
4.78(6) GeV from Ref. 0, evolve a s 5 \g) from g = M z 
to the matching scale g thr = £M t MC in the n/ = 5 ef¬ 
fective theory using coupled QCD x QED beta functions 
through four loops in QCD [7] and three loops in QED 
F221 , and evaluate there the MS couplings of the full SM 
from 


5 2 (M) 

g 2 (g)+g /2 (g) 

a(m) 

2//(m) 

gl{g) 


2 5 / 2 G F M 2 v [l + 6 w (g)], 
2 5 / 2 G F M z [l + 5 z (g)], 
2- 1 / 2 G F M 2 H [l + S H (g)}, 
2 3 / 4 G 1 J 2 Mf[l + S f (g)}, 
47rai 5) (^)[l + S aB {g)], (4) 


including the appropriate terms of 0(a n ) with n = 1,2 
[121 ill], O(oa s ) 0IMI1], and 0(a ") with n = 1, 2, 3,4 
[MS]- The threshold corrections 5i(g) in Eq. Q are 
expressed in terms of the MS couplings a(g) and a s (g), 
and the pole masses Mi. To avoid the theoretical un¬ 
certainty due to the hadronic contributions to the rela¬ 
tionship between a(g) and axh [2], we replace the lat¬ 
ter by G f in the set of basic parameters by extract¬ 
ing a(g) self-consistently from the exact relationship 
l/[47ra(/z)] = l/g 2 (g) + \/g' 2 {g) [H]. We stress that the 
MS couplings in Eq. Q are manifestly gauge indepen¬ 
dent and, thanks to partial tadpole cancellations, also 
finite in the limit Mh —> 0 [12] H4] . The pole masses 
M t and M b are subject to renormalon ambiguities of 
0(Aqcd)> which, for M t , are still small against the exper¬ 
imental error [2] and, for M b , are inconsequential because 
of the smallness of yb{g)- The use of MS masses m q (g) 
would avoid renormalon ambiguities at the expense of 
introducing unscreened tadpole contributions to restore 
gauge independence |9), which coincidentally reduce the 
scheme dependence of m t (g) [20, but spoil the pertur¬ 
bative expansion for m b {g) [12]. For completeness, we 
also study the MS mass parameter of the Higgs poten¬ 
tial, m 2 {g) = —2m%(g) = 2v 2 (g)X(g), using 

v 2 (g)=2- 1 / 2 Gr 1 [l + Af(g)], ( 5 ) 

where A r(g) to 0(a n ) with n = 1,2 and 0(aa s ) may 
be found in Ref. pH]. A f{g) is gauge independent, but 
diverges for Mh —t 0 due to unscreened tadpole contri¬ 
butions. We estimate the theoretical uncertainties in the 
MS parameters for £ = 1 due to unknown higher-order 
corrections by considering both scale variations and trun¬ 
cation errors. In the first case, we in turn put £ = 1/2 and 
2 in Eq. Q, return to £ = 1 using the RG equations, and 
select the larger one of the two deviations thus generated. 
In the second case, we find the full set of MS parameters 
for £ = 1, including besides those in Eq. Q also m, (g) 
with i = W, Z, H, f and v(g), by self-consistently solving 
the system of equations that express G F and M, entirely 
in terms of these parameters, so that unscreened tadpole 
contributions have to cancel numerically. We cast our 
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X 

x 0 

Ax as 

A xm w 

A x Mh 

A XMt 

Px 

<5a+ 

<^^tru 

9 

0.35838 

-3.8 x 10 e 

-2.3 x 10 -4 

-2.5 x 10~ b 

+7.1 x 10 t> 

+2.1 x 10 -3 

8.5 x 10" 3 

6.4 x 10~ 4 

g 1 

0.64812 +8.5 x 10~ 7 

+1.2 x 10“ 4 

-6.6 x 10~ 7 

-9.8 x 10 -6 

-5.2 x 10~ 3 

5.8 x 10" 5 

1.0 x 10” 3 

9s 

1.16540 +2.7 x 10" 3 

+8.9 x 10“ 8 

+7.8 x 10~ 8 

-4.0 x 10 -5 

-7.2 x 10 -2 

5.6 x 10" 5 


yt 

0.93517 

-3.6 x 10~ 4 

-1.3 x 10 -7 

-8.6 x 10~ 6 

+5.1 x 10 -3 

-5.2 x 10" 2 

8.0 x 10 -4 

1.2 x 10" 3 

Vb 

0.01706 

-5.7 x 10" 5 

-5.1 x 10 _1 ° 

+1.3 x 10~ 7 

-2.4 x 10 -7 

-9.2 x 10 -4 

2.5 x 10" 4 

1.1 x 10 -3 

X 

0.12714 

-6.2 x 10“ 6 

-4.2 x 10 -7 

+8.2 x 10~ 4 

+6.4 x 10 -5 

-2.0 x 10 -2 

5.8 x 10" 4 

5.5 x 10” 4 

m 

131.86 

-2.6 x 10 -3 

-4.4 x 10 -4 

+3.8 x lO^ 1 

+1.2 x 10 1 

+2.6 

7.3 x 10" 1 

4.1 x 10 -2 


TABLE I: Coefficients in Eq. |6|. The entries in the last row are given in units of GeV. 


results for x = g, g',g a ,yt, Vb, A, m in the form 

af\M z ) ai 5) ’ exp (M z ) 


v(g) = x 0 + Ax a 


A ai 5hexp (M z ) 


M w -M™ p , A M H -M e H xp 

+ ^M W AM exp + A XMh AM exp 


L W 

M t - MT P 


,thr 


i a lv±t i o A 4 A 4 ir 

+ AX Mt --h Px -Iu7— ± 0X L 


AM t exp 


,thr 


± Sx t 


( 6 ) 


allowing for linear extrapolations in the least precisely 
known input parameters quoted above [2] and /i thr , which 
we disentangle from M t exp = M t MC , and list the coeffi¬ 
cients in Tabic |U 

We now in turn apply criterion © and the approach 
of Ref. jJO] and write the resulting critical masses and 
associated scales X = M+ 1 , gf' . M'? rl , gf' with i = t,H 
in the form 


X = Xn + AX n 


af\M z ) a^jMz) 
A ai S) ’ exp (M z ) 


M - M exp 

' AX m — XT7 - ± 6X , 

AM exp 


pai + SX^±6X tIU , (7) 


where M = Mp (M t ) if i = t ( H ), AA' ai , and AXm are 
the lcr errors due to ai 5 ' 1 (M z ) and M, respectively, <LY par 
are the residual parametric errors combined in quadra¬ 
ture, 5X± are the shifts due to the choices £ = 2 ±1 , and 
<5X tru are the truncation errors induced by those in Ta¬ 
ble [I] The coefficients in Eq. |7]) are collected in Table [TT| 
M t crl is 0.20 GeV larger than M t crl , and Mp 1 is 0.40 GeV 
smaller than M p'. These shifts reflect the scheme depen¬ 
dence. gf' and g c p fall slightly short of Mp, for which 
log w Mp = 19.086, where the SM definitely ceases to 
be valid, while gf' and gff lie appreciably beyond M Pl 
which is an inherent problem of Ref. [20 ] and was cured 
there by the ad hoc introduction of some new dimension- 
six operator. In the remainder of this Letter, we concen¬ 
trate on the approach based on Eq. © 0- 

To assess the significance of the higher-order correc¬ 
tions that were not yet included in Ref. [5], namely the 
full 0(a 2 ) terms in 5fg) with i = W,Z,H,q ]2. Til , 
the 0(aa s ) term in S as (g) T3], and the 0(a 4 ) terms 
in S as {g) ijj! and 8 q {g) PSj, we switch them off one at 


a time. The resulting central values and scale depen¬ 
dencies of the critical parameters are also contained in 
Table |© The 0(a 2 ) terms in dfg) J2'„ 03] shift Mf' 
and Mp' by —0.11 and +0.24 GeV, respectively, and re¬ 
duce their scale uncertainties by almost a factor of 3. 
On the other hand, the 0(a 4 ) terms in 5 q (g) [18] pro¬ 
duce larger and opposite shifts in Mf' and Alp', namely 
+0.20 and —0.42 GeV, respectively, but merely reduce 
their scale uncertainties by less than 10%. The 0(aa s ) 
0 and 0(a 4 ) [T7] terms in 5 as (g) are much less signifi¬ 
cant. All these observations approximately carry over to 
Mf' and Mf. 

Apart from the issue of gauge dependence, our analysis 
differs from that of Refs. mm in the following respects. 
In Refs, mm, the 0(aa s ) term in 6 a „{g) [12 and the 
0(a 4 ) terms in 8 as (g) [T7] and 8 q (g) [IS] were not in¬ 
cluded; g thr was affected by the M t MC variation, which 
explains the sign difference in the corresponding shift in 
Mp'; and the scale uncertainties were found to be ap¬ 
proximately half as large as here for reasons unknown to 
us. 

In Fig. [TJ the RG evolution flow from g thr to g crl and 
beyond is shown in the (A,/?a) plane. The propagation 
with g of the lcr and 3cr confidence ellipses with respect 
to ALA 10 and Mp tells us that the second condition in 
Eq. © is almost automatic, the ellipses for g = 10 18 GeV 
being approximately degenerated to horizontal lines. For 
default input values, A (g) crosses zero at g = 1.55 x 
10 10 GeV. The contour of M f approximately coincides 
with the right envelope of the 2<r ellipses, while the one 
of Mp', which relies on Mf IC , is driven outside the 3<r 
band as g runs from gf to g thT . 

Our upgraded and updated version of the familiar 
phase diagram 1® HI EDI Ml is presented in Fig. [© 
Besides the boundary of the stable phase defined by 
Eq. ©, on which the critical points with Mf and 
Mp' are located, we also show contours of A (g°) = 0 
and f3\(g°) = 0. The demarcation line between the 
metastable phase and the instable one, in which the life¬ 
time of our vacuum is shorter than the age of the Uni¬ 
verse, is evaluated as in Ref. [2D] and represents the only 
gauge-dependent detail in Fig. [2l The customary confi¬ 
dence ellipses with respect to Mf AC and Mp, which are 
included Fig. ©for reference, have to be taken with cau¬ 
tion because they misleadingly suggest that the tree-level 
mass parameter M t MC and its error [2] identically carry 
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X 

A'o 

ax Q6 

A X M 

dApar 

SX+ 

sx~ 

sx tIU 

s° (c,2) 

aO(ctCt s , 0 ^ ) 
Ga s 

rO{a*) 

O q 

Mr 

171.44 

0.23 

0.20 

0.001 

-0.36 

0.17 

-0.02 

171.55 ;?;?? 

171.43;°;?? 

171.24;?;?? 

log 10 Mt rl 

17.752 

-0.051 

0.083 

0.007 

0.007 

-0.006 

-0.002 

1 7 7QQ+0.062 

1 l . 1 oo _ 0 008 

i 7 7 ^ 4 + 0.007 
i /. /O4_o 006 

1 7 7^1 +0-UU7 
l (. 1 dl_ 0 007 

M/P 

129.30 

-0.49 

1.79 

0.002 

0.72 

-0.33 

0.04 

129.061??? 

129 32 +0,73 

129.721?;?? 

logm^S 1 

18.512 

-0.158 

0.381 

0.008 

0.173 

-0.082 

0.008 

18.4951?;??? 

18.518!?;??? 

18.6021?;??? 

Mr 

171.64 

0.23 

0.20 

0.001 

-0.36 

0.17 

-0.02 

171-74;?;?? 

171.63;°;?? 

171.43;?;?? 

logmAf 1 

21.442 

-0.059 

0.094 

0.005 

-0.083 

0.022 

0.002 

21.485; 0 0 ;?? 3 5 

21.445;°;??? 

91 441 -0.072 

ZI.44:l +0 Q14 

M/P 

128.90 

-0.49 

1.79 

0.003 

0.73 

-0.34 

0.04 

128.671??? 

128.921??? 

129.321?;?? 

logio AS 1 

22.209 

-0.181 

0.436 

0.007 

0.092 

-0.062 

0.013 

22.2011°-“® 

99 91 7+0-094 
< -0.062 

99 *31 O+0.113 
082 


TABLE II: Coefficients in Eq. (M) and central values with scale dependencies obtained upon switching off the 0(a 2 ) terms in 
with i = W, Z, H, q, the 0[aa a ) and O(o?) terms in 5 as (fJ-), and the 0(a a ) terms in S q (fj.) one at a time. The unit of 
mass is taken to be GeV. 



FIG. 1: RG evolution of X(y) from /i thr to /i crl and beyond 
in the (A, / 3 a ) plane for default input values and matching 
scale (red solid line), effects of la (brown solid lines) and 3 a 
(blue solid lines) variation in M t MC , theoretical uncertainty 
due to the variation of £ from 1/2 to 2 (upper and lower 
black dashed lines with asterisks in the insets), and results 
for M t crl (green dashed line) and M/p (purple dashed line). 
The la (brown ellipses) and 3cr (blue ellipses) contours due to 
the errors in M t MC and Mh are indicated for selected values 
of y. The insets in the upper right and lower left corners refer 
to fj, = A7 f MC and /r = 1.55 x 10 10 GeV, respectively. 


over to M t , which is actually the real part of the complex 
pole position upon mass renormalization in the on-shell 
scheme (251 . In view of the resonance property, a shift of 
order T t = 2.00 GeV [2] would be plausible, which should 
serve as a useful error estimate for the time being. 

In conclusion, we performed a high-precision analy¬ 
sis of the vacuum stability in the SM incorporating full 
two-loop threshold corrections [51 EMI], three-loop beta 
functions jl], and 0(a?) corrections to the matching and 
running of g s [3 E3 and y q 5J [15], and adopting two 
gauge-independent approaches, one based on the criti¬ 
cality criterion p]) for A(/r) 5] and one on a reorgani¬ 
zation of V e s(H) so that its minimum is gauge inde¬ 
pendent order by order [20] . For the M t upper bound 



FIG. 2: Phase diagram of vacuum stability (light-green 
shaded area), metastability, and instability (pink shaded area) 
in the ( Mh , M t ) plane, contours of A(/r°) = 0 for selected val¬ 
ues of fjP (purple dotted lines), contours of = 0 for se¬ 

lected values of (solid parabolalike lines) with uncertainties 
due to la error in a+(Mz) (dashed and dot-dashed lines), 
critical line of Eq. (solid green line) with uncertainty due 
to la error in cri+Mz) (orange shaded band), and critical 
points with M t crl (lower red bullet) and M/p (right red bul¬ 
let). The present world average of (M t MC ,Mrr) (upper left 
red bullet) and its lcr (purple ellipse), 2 a (brown ellipse), and 
3 <j (blue ellipse) contours are marked for reference. 


we thus obtained M t crl = (171.44 ± 0.30p°gg) GeV and 

M f crl = (171.64±0.30^o'gg) GeV, respectively, where the 
first errors are experimental, due the lcr variations in the 
input parameters [2], and the second ones are theoretical, 
due to the scale and truncation uncertainties. In want of 
more specific information, we assume the individual error 
sources to be independent and combine them quadrati- 
cally to be on the conservative side. The 0.20 GeV dif¬ 
ference between the central values of M/ rl and M t crl in¬ 
dicates the scheme dependence, which arguably comes 
as a third independent source of theoretical uncertainty. 
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As our final result, we hence quote the combined value 
M t cri = (171.54 ± 0.30 ^Q 4 ^) GeV, which is compatible 
with = (173.21 ± 0.87) GeV at the 1.3<r level. In 
view of this and the present lack of knowledge of the pre¬ 
cise relationship between and M t MC and M t mentioned 
above, the familiar notion Emm that our vacuum is 
metastable is likely to be premature [24]. 
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DFG Grant No. SFB 676, MES of Russia Grant No. MK- 
1001.2014.2, the Heisenberg-Landau Programme, and 
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